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We present a detailed phenomenological analysis of single-inclusive hadron production at the 
CERN-LHC in both proton-proton and proton-lead collisions. First data from the LHC exper- 
iments on charged hadron spectra are compared to next-to-leading order QCD expectations, and 
predictions are made for identified pion, kaon, and proton distributions differential in transverse mo- 
mentum and rapidity for LHC energies from 900 GeV to f 4 TeV. The results are obtained with the 
latest sets of vacuum fragmentation functions based on global QCD analyses, and recently proposed 
medium modified fragmentation functions are used to model hadronization in proton-lead collisions 
assuming standard QCD factorization. Besides estimating theoretical ambiguities due to the choice 
of factorization and renormalization scales and parton densities, we carefully assess uncertainties due 
to our present knowledge of parton-to-hadron fragmentation functions with the Lagrange multiplier 
technique, ft is outlined to what extent future LHC data will contribute to further our quantitative 
understanding of hadronization processes. 
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I. INTRODUCTION AND MOTIVATION 

Single-inclusive hadron production is increasingly be- 
coming a precise and versatile tool to quantitatively 
study various aspects of Quantum Chromodynamics 
(QCD), supplementing the insights already gained from 
fully inclusive observables. First and foremost, precise 
data obtained in electron-positron, lepton-nucleon, and 
hadron-hadron collisions at different center-of-mass sys- 
tem (c.m.s.) energies yS have established collinear fac- 
torization []J as the foundation for a perturbative QCD 
(pQCD) approach to hard scattering processes involving 
identified hadrons produced at large transverse momen- 
tum pt ■ 

Within the framework of collinear factorization, the 
presence of a hard scale pt allows one to compute with 
remarkable phenomenological success various kinematic 
distributions of the observed final-state hadron H as a 
convolution of partonic short-distance scattering cross 
sections, calculable as a perturbative expansion in the 
strong coupling a s , and universal but non-perturbative 
functions capturing the long-distance physics as repre- 
sented by parton distribution and fragmentation func- 
tions. 

By now, information on these non-perturbative inputs 
is routinely extracted in "global QCD analyses" by con- 
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sistcntly comparing data sets from many different exper- 
iments and processes with theoretical expectations at a 
given order in perturbation theory. In particular for par- 
ton density functions (PDFs) , such type of analyses have 
reached a high level of sophistication in recent years @- 
4]. Various methods have been put forward to arrive 
at faithful estimates of the remaining uncertainties in 
our understanding of the hadronic structure and 
how they affect, e.g., the level at which we have control 
of standard model background processes at the CERN 
Large Hadron Collider (LHC). 

Extractions of parton-to-hadron fragmentation func- 
tions (FFs) and their uncertainties are a much more re- 
cent achievement 043 > mainly hampered by the con- 
siderably smaller amount of experimental results avail- 
able. In general, reducing uncertainties of FFs is of 
utmost importance to further our understanding of the 
hadronization process, test the scale dependence of FFs 
as predicted in pQCD, and to delineate kinematic regions 
where factorized, color- independent, collinear parton-to- 
hadron fragmentation provides a sensible approximation. 

The limited knowledge of FFs has a serious impact on 
different studies, such as, for instance, global QCD anal- 
yses of data taken in collisions of longitudinally polarized 
nucleons and leptons aiming to address the fundamental 
question of how the spin of the nucleon is composed of 
the intrinsic spins and orbital angular momenta of quarks 
and gluons [lfjj. Since many of the data available so 
far involve identified hadrons in the final-state, either in 
semi-inclusive deep-inelastic lepton-nucleon or in proton- 
proton scattering processes, FFs are a crucial ingredient 
of theoretical analyses of the spin structure of the nu- 
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cleon. Processes with identified hadrons are also increas- 
ingly used as decisive probes to gain insight about the 
properties and nature of cold and hot nuclear matter in 
heavy ion collisions both at the BNL Relativistic Heavy 
Ion Collider (RHIC) and soon also at the LHC. Again, 
precise knowledge of FFs and possible modifications in- 
duced by the nuclear medium are crucial ingredients that 
have received increasing attention recently 

In the following, we will assess in some detail how up- 
coming data from the LHC on high hadron production 
in an unprecedented kinematic regime will help to further 
our knowledge of FFs both in the vacuum and in a nuclear 
environment. To this end, we first need to discuss briefly 
some theoretical preliminaries such as the formalism for 
single-inclusive hadron production within pQCD factor- 
ization, the range of applicability of universal vacuum 
FFs, estimates and propagation of uncertainties, and a 
recent proposal for medium modified FFs. 

A. Framework for hadron production at the LHC 

Throughout this paper, we are interested in the single- 
inclusive invariant cross section for the production of a 
hadron H with energy E and momentum p in hadron- 
hadron or hadron-nucleus collisions. Assuming pQCD 
factorization, the relevant theoretical expression for pp 
collisions at the LHC schematically reads 

d 3 a H 

a.b,c 

® da ab ^ cX (S, a s ,x a ,x b ,z c ,fif,iJ, r , {if >),(!) 

where the sum is over all contributing partonic subpro- 
cesses ab — > cX contained in the perturbatively calcula- 
ble short-distance scattering cross sections da a b^cX- 

The scales fif and fif are introduced to factorize ini- 
tial and final-state collinear singularities into the scale 
dependent PDFs and FFs, f a ,b{x a ,b, M/) an d (z c , (J,f>), 
respectively. \i r denotes the energy scale at which a s is 
being renormalized. The residual dependence of Eq. ([T]) 
on the arbitrary scales |tt/,/',r can be taken as a measure 
for the theoretical ambiguity due to the truncation of 
the perturbative series at a given fixed order in a s . To 
estimate its impact, we will follow the usual procedure 
and vary the scales within a factor of two around the de- 
fault choice pt- It turns out that scale variations are the 
dominant theoretical uncertainty for hadron production 
at the LHC. 

x a ^ are the fractions of longitudinal momentum of the 
colliding hadrons taken by the interacting partons a and 
b. Similarly, z = z c denotes the collinear momentum frac- 
tion of the fragmenting parton c carried by the produced 
hadron H. Neither x a ^ nor z are measurable quanti- 
ties, and any given data point characterized by the px 
of the hadron H and the c.m.s. energy ^J~S probes both 
the PDFs and the FFs at a different range of momentum 
fractions and scales of 0(px)- Likewise, different hadron 



species H = 7r° , , K , . . . determine complementary 
aspects of the hadronization process, i.e., different FFs 
. The relevance of each parton flavor c depends on 
the quark content of H and on px and which con- 
trol the contributions of the various partonic channels 
ab —> cX to the sum in Eq. (J]). We will demonstrate 
how the different partonic subprocesses arc expected to 
contribute to hadron production yields at LHC energies. 

A noteworthy property of the invariant cross section 
in Eq. ([1} is its approximate power-law behavior [l2l . Il3j 

a inv = E^-=F(x T )/p T (xT '^ (2) 

for fixed xt = %Pt /VS '. In the naive, scale-invariant par- 
ton model one expects scaling with n = 4. The running 
of a s and the scaling violations of the PDFs and FFs as 
predicted by pQCD lead to deviations from exact scaling 
for Eqs. JT]) and @, i.e., n = ti(xt,V^S), which can be 
explored by comparing, e.g., xt hadron spectra at differ- 
ent \/S. We shall briefly touch upon i T scaling at LHC 
energies and estimate theoretical scale uncertainties in 
predicting n = u(xt, V~S)- 

B. Applicability of FFs and uncertainty estimates 

Since the LHC sets a new energy frontier for hadron 
production, it is crucial to first convince ourselves that 
the accessible range of z is still compatible with the ap- 
plicability of the concept of FFs [14J within the factor- 
ized framework as outlined above. Contrary to PDFs, 
where the small x a .b regime is rather well explored down 
to momentum fractions of 0(1O~ 4 ) [2 4] and has been 
subjected to a very detailed theoretical scrutiny, the phe- 
nomenological access to FFs is much more restricted. Ne- 
glected hadron mass effects, potential higher twist cor- 
rections, and instabilities in the timelike scale evolution 
limit the usage of FFs to rather large values of z > 0.05; 
see, e.g. Ref. Q- We shall show that for all practical 
applications at the LHC the condition z > 0.05 is well 
met as the bulk of the cross section for inclusive hadron 
production samples on average large momentum frac- 
tions, (z) ~ 0.5, well above the kinematic lower limit 
z w 2 PT /VS. 

Upcoming experimental results for single-inclusive 
hadron production at the LHC can be straightforwardly 
included [lj| in existing global QCD analyses of FFs 
0, [H , which are the most efficient and consistent method 
to deconvolute information on these non-perturbative 
functions from the interplay of various observables mea- 
sured at different energy scales. The large range of trans- 
verse momenta pt accessible at the LHC will allow for un- 
precedented studies of evolution effects for FFs. We note 
that for both the timelike scale evolution of FFs [l6| and 
the partonic hard scattering cross sections da a b^ c x [HI 
in Eq. ([T]), pQCD calculations at next-to-leading order 
(NLO) accuracy are "state-of-the-art" and used through- 
out this work. In any case, they are mandatory for an 
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accurate and meaningful comparison of theory and data 
due to often sizable NLO QCD corrections and, in par- 
ticular, to allow for estimates of scale uncertainties. 

As mentioned above, besides establishing a small set 
of "best-fit" parameters in a x 2 minimization to model 
the functional form of the FFs for different flavors and 
hadrons, assessing their uncertainties and propagating 
them to physical observables is an equally important 
goal. Here, the most robust technique is based on La- 
grange multipliers 043 which makes no assumptions 
about the behavior of the x 2 profile near its minimum. 
We illustrate the usefulness of this method by estimating 
the uncertainties from FFs for charged hadron produc- 
tion at the LHC and Fermilab's Tevatronpp collider. For 
the latter, data from the CDF collaboration [18j have re- 
cently caused some stir as the measured cross section 
above pr > 20 GeV exceeds theoretical expectations by 
orders of magnitude Since the data are also in ex- 
cess of single-inclusive jet cross section measurements, 
the most likely explanation seems to be some experimen- 
tal problem [ijj. Nevertheless, we believe it is an useful 
exercise to evaluate the theoretical ambiguities caused by 
FFs also in this case. It turns out that naive estimates, 
for instance, by comparing the results obtained with two 
different sets of FFs, often seriously underestimate un- 
certainties from FFs. 



C. FFs in a nuclear environment at the LHC 

It is well known that results for hadron production 
processes occurring in a nuclear medium can differ sig- 
nificantly from similar experiments involving only light 
nuclei or proton targets, showing both suppression and 
enhancement of the rates depending on the details of the 
observable. Recent examples include production rates 
of pions and kaons in semi-inclusive deep-inelastic scat- 
tering off different nuclei as provided by the HERMES 
experiment |20| and in deuteron-gold collisions measured 
at RHIC [21|. The origin of the observed nuclear mod- 
ifications has been attributed to a variety of conceiv- 
able mechanisms and models [ll[ besides the well-known 
modification of parton densities in nuclei (nPDFs) [22i — 
l24j |. Available models incorporate ideas based on inter- 
actions between the nuclear medium and, e.g., the final- 
state hadron or the seed partons before the hadronization 
takes place, and reproduce, with different degree of suc- 
cess, some features of the data; for recent reviews, see 
Ref. [HI- 

nPDFs provide an effective and phenomenologically 
successful way to factorize the influence of the nuclear 
environment on the interacting partons into sets of uni- 
versal functions which scale in energy like ordinary PDFs 
and can be obtained in global QCD fits to available data 
[22U2~i ] . The quite natural extension of this idea to final- 
state nuclear effects has been put forward only very re- 
cently by introducing the concept of medium modified 
fragmentation functions (nFFs) First QCD fits for 



identified pions and kaons were provided in Ref. [25j re- 
cently. As for nPDFs, it was demonstrated that within 
the precision of the available data, universal nFFs are 
a viable concept and factorization similar to Eq. ([T]) 
holds at least approximately [25[ despite being much 
more speculative than in pp collisions (26j . The nuclear 
A dependence of the nFFs can be most economically 
parametrized in a convolution approach [25j 

D? /A (z,no)=J ^-W?(y,A,(, )D? (j,/i ) (3) 

which relates the nFFs D^ A to the fairly well known 
vacuum FFs of DSS [7] at some initial scale through 
a weight function Wjr with only a small amount of extra 
parameters. nFFs at scales \x > /xo are then obtained by 
applying the standard timelike evolution equations [l6| . 

Combined with nPDFs, nFFs allow one to treat a 
large class of hard hadron production processes where 
a nucleus collides with a lepton or a nucleon (light nu- 
cleus) in a consistent pQCD framework based on factor- 
ization. Exploiting the predictive power of the factorized 
approach, we will provide predictions for pion production 
in future proton-lead (pPb) collisions at the LHC in a 
wide range of pr and rapidity. pPb collisions at the LHC 
are conceivable up to a c.m.s. energy of VS ~ 8.8 TeV 
[27l ] though not part of the initial LHC physics program. 
In addition, we shall discuss how the admixture of the dif- 
ferent contributing partonic subprocesses ab — ¥ cX and 
the range of probed momentum fractions z is expected to 
change in a nuclear environment at LHC energies. Such 
measurements will be crucial to further our knowledge 
of hadronization in a nuclear medium by exploring to 
what extent factorization breaking effects due to inter- 
actions of partons with the medium come into play and 
limit the usefulness of nPDFs and nFFs. Comparisons 
of hadron rates obtained pp and pPb collisions will also 
help to unravel and understand the properties of hot and 
dense QCD matter. 

The remainder of the paper is organized as follows: in 
the next Section, we compare first results for px spectra 
of unidentified charged hadrons from ATLAS [28| and 
CMS (2t| with pQCD calculations at NLO accuracy. In 
Sec. Ill, we present expectations for both identified and 
unidentified hadron production cross sections in pp col- 
lisions at the LHC in a broad range of pt and rapidity 
y. The results are supplemented by studies of theoreti- 
cal scale and PDF ambiguities, the relevance of different 
partonic subprocesses, and the ranges of momentum frac- 
tions Xa t b and z predominantly probed. We touch upon 
xt scaling and give estimates of FF uncertainties for 
charged particle yields at the LHC and the Tevatron us- 
ing the Lagrange multiplier technique. Pion production 
in pPb collisions and medium modified FFs are discussed 
in Sec. IV. We summarize our main results in Sec. V. 
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II. COMPARISON TO FIRST LHC DATA 

Before turning to a detailed discussion of theoretical 
expectations for single-inclusive hadron production at 
LHC energies in Sec. JIII1 we take a brief look at first 
results from ATLAS [2^ and CMS [H| for p T differen- 
tial charged hadron yields. 

Charged hadron multiplicities are the first results 
of the LHC physics program, and data were reported 
shortly after the startup of the LHC by the ALICE 
[ISIlil, ATLAS [H], and CMS ^ experiments which 
impressively demonstrated the readiness of their detector 
systems. They provide measurements of the number of 
charged hadrons N c h with respect to (w.r.t.) their trans- 
verse momentum px and pseudorapidity i[. The main 
result so far is the increase of the pseudorapidity density 
dN c h/dr] at central rapidities with the c.m.s. energy of 
the pp collisions in the range \/S = 0.9 - 7 TeV 0, |30[. 

The pseudorapidity density dN c h/drj is not amenable 
to pQCD calculations based on factorization, Eq. ([TJ, 
as the bulk of the produced hadrons has very low pt, 
well below 1 GeV, and a hard scale is lacking. However, 
both ATLAS @ and CMS [H| also present data for 
dN c h/dpT in the rapidity range |?7| < 2.5 and \r]\ < 2.4, 
respectively, at sufficiently high pt- Very recently, also 
ALICE published a pt spectrum at ^/S = 900 GeV and 
\v\ < 0.8 [3l|, but data tables are not yet available for 
comparison to theory. In each case, the charged parti- 
cle multiplicities are normalized to the number of inelas- 
tic non-single-diffractive (NSD) interactions but based on 
slightly different event selections. 

In order to compare the obtained invariant yields 
\/(2i:pT)d 2 N c h/d'qdpT with pQCD predictions for 
Ed 3 a H /dpm Eq. (TTJ), one needs to normalize them with 
the total NSD cross section CAf.^n , which unfortunately is 
not specified in Refs. 0, [2||[3lj]. However, we estimate 
vnsd for each available data set by scaling the theo- 
retical single-inclusive hadron yields at NLO by 1/(Jnsd 
and fitting a common, i.e., p^-independent, value for each 
experiment. We use the CTEQ6.6 set of PDFs [2j and 
the DSS parton-to-unidentified charged hadron FFs [8j in 
Eq. (fTJ) and choose [if = [if = Mr = A* = Pt- The results 
of the NLO calculations are shown as solid lines in Figs. [I] 
and EJand compared to ATLAS [H[ and CMS [HI data, 
respectively. Note that only data with px > 1 GeV are 
displayed. We refrain from giving theoretical expecta- 
tions based on Eq. (TTJ for pt < 1.5 GeV where power 
suppressed corrections to the factorized pQCD frame- 
work and other non-perturbative soft contributions are 
relevant. In any case, all sets of PDFs Q and FFs Q are 
not applicable for too small scales [i ~ pr- 

The shaded bands in Figs. [T] and [5] give an indication 
of the theoretical uncertainties due to the truncation of 
the perturbative series in Eq. (TJJ) at NLO accuracy. As is 
customary, they are obtained by simultaneously varying 
all scales in Eq. ([T]) by a factor of 2 around the default 
choice fi r = pt- Given the fact that the normalization of 
the theoretical results has to be determined by a fit, we 




FIG. 1: Comparison of the single-inclusive charged hadron 
yield per NSD event in the rapidity range 77 < 2.5 in pp 
collisions at y/S = 900 GeV measured by ATLAS [13] with a 
NLO pQCD calculation using the DSS FFs [| and CTEQ6.6 
PDFs 2] (solid line) . The shaded band indicates the theoret- 
ical uncertainty from varying the scales in Eq. (ft) by a factor 
of 2 around the default choice [if = [iti — ji r = pr- Note 
that the normalization of the curve is determined by a fit; see 
text. 



refrain from studying other, usually subleading, sources 
of uncertainties like variations of PDF sets, at this point. 

The overall agreement between the px dependence as 
predicted by pQCD calculations at NLO accuracy and 
data is excellent, and the drop of the cross sections with 
increasing px over many orders of magnitude is accu- 
rately reproduced. As was already noticed in pp collisions 
at BNL-RHIC at lower c.m.s. energies, NLO pQCD cal- 
culations at the leading-twist level tend to describe sin- 
gle inclusive hadron yields well down to fairly low values 
of pt — 1.5 GeV. We note a mild tension with CMS 
data [H in Fig. CD for y/S = 7 TeV and p T > 4.5 GeV, 
which are all at the lower edge of the scale uncertainty 
band. Notice that of all hadron production data, the 
ones taken at V~S — 7 TeV probe the smallest x a ^b values 
in the PDFs so far, down to a kinematic limit of about 
Xt — 5 • 10~ 4 . 

Upcoming precision data from the LHC will decisively 
test all theoretical aspects of inclusive hadron production 
in an unprecedented energy range, from non-perturbative 
PDFs and FFs to the validity of the leading-twist approx- 
imation in Eq. (TTJ. In addition to unidentified charged 
hadron yields measured so far, the ALICE experiment 
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FIG. 2: As in Fig. [T] but now in the rapidity range < 2.4 
and for three different c.m.s. energies = 0.9, 2.36, and 
7TeV as measured by the CMS experiment [29 ]. 



is designed to identify various hadron species like pions, 
kaons, and protons at central rapidities. These results 
will provide a vital input for future global analyses of 
FFs. Detailed theoretical expectations for single inclu- 
sive hadron production at the LHC will be discussed in 
the next Section along with estimates of uncertainties. 



III. EXPECTATIONS FOR SINGLE-INCLUSIVE 
HIGH-pt HADRON PRODUCTION IN pp 
COLLISIONS AT TeV ENERGIES 

Encouraged by the first results from the LHC experi- 
ments and the successful comparison to theoretical expec- 
tations based on the DSS set of fragmentation functions 
@, 111 , we now turn to a detailed discussion of both iden- 
tified and unidentified hadron production cross sections 
to be studied at the LHC in the future. While ATLAS, 
CMS, and ALICE can all measure unidentified charged 
hadron yields, only ALICE is capable of tagging different 
hadron species at central rapidities, which is crucial for 
detailed tests of our current knowledge of FFs and the 
underlying framework of pQCD. 

Our results are supplemented by estimates of theo- 
retical scale, PDF, and FF uncertainties. To elucidate 
the impact of future LHC data on our understanding of 
pQCD factorization and, in particular, hadronization, we 
discuss the relevance of different partonic subprocesses 
dtfab^cX in Eq. (p}, the role of quark and gluon frag- 
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FIG. 3: Upper panel: NLO cross section for tt + production at 
four different c.m.s. energies, integrated in the pseudorapidity 
range \r]\ < 1.5. The shaded bands indicate the uncertainty 
associated with the variation of the scales in the range pr/2 < 
f^fj'.r < 2pr- Lower panel: ratio of the ir~ and n + yields. 



mentation, and the ranges of momentum fractions which 
will be predominantly probed. 

If not stated otherwise, we use the CTEQ6.6 set of 
PDFs and the associated Hessian sets for PDF uncer- 
tainty estimates @ • The parton-to-hadron FFs are taken 
from DSS 0, H|, and uncertainties related to FFs will 
be estimated with the Lagrange multiplier method 0- 
8] . All calculations are performed at NLO accuracy with 
hard scattering cross sections, PDFs, and FFs in Eq. (Q} 
evaluated in the MS scheme. We take the transverse mo- 
mentum pt of the produced hadron as the default choice 
for the factorization and renormalization scales, /**/./' and 
/i r , respectively, in Eq. ([TJ. To estimate theoretical ambi- 
guities associated with the truncation of the perturbativc 
series at NLO accuracy we vary all scales by a factor of 
two up and down their central value px as is commonly 
done. 



A. Pion, kaon, and proton production at the LHC 

We start off with showing the energy dependence of the 
Pt differential cross section for charged pion production 
integrated in the pseudorapidity range I77I < 1.5 in the 
upper panel of Fig. for y/S = 0.9, 2.36, 7, and 14 TeV. 
The shaded bands give an estimate of the theoretical 
ambiguities due the variations of fJ-fj>, r as described 
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FIG. 4: Upper panel: comparison of the yields for tt + , 
K + , and p at NLO accuracy for pseudorapidity rj = and 
i/S — 7TeV. As before, the shaded bands give an estimate 
of the theoretical scale ambiguity. Lower panel: PDF un- 
certainty for 7r + production w.r.t. the solid line in the up- 
per panel, estimated with the Hessian eigenvector sets for 
CTEQ6.6 @] (solid band) and compared to the ambiguity 
from scale variations (hatched band). 



above. The typical scale uncertainty for \/~S — 14TeV 
amounts to about 25% at px ~ 20 GeV and 14% at 
Pt ~ 100 GeV. At the same values of px but lower 
c.m.s. energies, = 0.9 TeV, the corresponding un- 
certainties can reach up to 30%. 

The lower panel shows the expected difference of tt + 
and 7r~ yields as predicted by the DSS set of FFs Q. 
We recall that at present, the charge separation of FFs is 
mainly constrained by semi-inclusive deep-inelastic scat- 
tering (SIDIS) data at a relatively low scale Q 2 ~ 2.5GeV 
and to a much lesser extent by data from RHIC Q. 
As a consequence, uncertainties are still sizable in cur- 
rent sets of FFs. As can be seen from Fig. [3J the ratio 
da'* I da'* drops only very slowly with increasing pt for 
both \fS = 14 and 7 TeV. This is due to the dominance 
of gluon initiated hard scattering and hadronization pro- 
cesses, which prevails up to fairly highp T ; see Figs. [5] (a) 
and (c) and the discussions below. Clearly, a precise mea- 
surement of da* I da* is of great phenomenological im- 
portance but at the same time also very challenging at 
nominal LHC energy as sufficiently large integrated lumi- 
nosities are required to resolve effects of a few percent at 
high px- Lower c.m.s. energies, = 2.36 and 0.9 TeV 
are more favorable as they probe larger momentum frac- 



tions x a fi in the nucleon in the same range of px, i.e., 
larger xt = 2px/y/S, where quark initiated scattering 
processes more and more dominate. 

The expectation that dcr 7r+ > da~* can be understood 
by inspecting the role of different partonic subprocesses 
in pp collisions; see Fig. |6] (a) below. Quark-gluon scat- 
tering is the second most important channel for hadron 
production. The abundance of u quarks in a proton and 
the ud valence flavor structure of a 7r + , i.e., > D„ , 
explains the observed hierarchy for the production 
yields found in Fig. [3j 

Figure [4] compares the cross sections for ir + , K + , and 
proton production as a function of px for ^/S = 7 TeV 
and 77 = 0. At pt ~ 20 GeV pions yields are roughly a 
factor of 3 larger than those for K + , and the production 
of protons is suppressed even further. Together, pions, 
kaons, and protons account for almost the entire yield 
of charged hadrons; see Fig. [8] The mixture of hadron 
species is largely independent of px, with the fraction of 
pions decreasing slightly as px increases. This is because 
the mean value of z changes only very slowly with px, 
see Fig. [6] (e). At higher z, fragmentation into heavier 
hadrons like kaons and protons is somewhat enhanced 
w.r.t. those into pions @,Hj]- 

The typical PDF uncertainty for hadron production at 
the LHC is illustrated in the lower panel of Fig. [4] and 
amounts to about 5% at px — 20 GeV. It is computed 
with the help of the 44 Hessian eigenvector sets provided 
by CTEQ [2| and compared to variations of the cross sec- 
tion due to the choice of scales Hfj>, r in Eq. (pj. Because 
of the dominance of gluon initiated scattering processes, 
see Fig. [5] (a), the PDF uncertainty reflects to a large 
extent the present ambiguity in the gluon PDF at the 
relevant scale fj, ~ px and range of momentum fraction 
x. In particular at smaller values of px, the theoretical 
ambiguity due to the truncation of the pQCD series at 
NLO turns out to be by far the most relevant one. 

The pseudorapidity dependence of the production 
cross section for neutral pions at \/~S = 7 TeV is shown 
in Fig. [5] for various fixed values ofpx- Again, the typi- 
cal scale uncertainty is indicated by the shaded bands for 
px = 10 and 50 GeV. Measurements of the 77 dependence 
for fixed px in a wide range are of great phenomenologi- 
cal interest as they emphasize different partonic subpro- 
cesses and momentum fractions x a ^ and z as compared 
to the px differential yields shown in Figs. |3]and|U 

These aspects are explained in some detail in Fig. [5] 
The panels show the relevance of the different NLO par- 
tonic subprocesses da a b (upper row), the relative frac- 
tions of quark and gluon fragmentation (middle row), 
and the mean values of the momentum fractions x a ,b and 
z (lower row) probed in tt° production at yS = 7 TeV. 
The results on the left hand side (l.h.s.) are given as a 
function of px for rj = 0, and the panels on right hand 
side (r.h.s.) are differential in r\ for fixed px = 25 GeV. 
For different c.m.s. energies, the plots roughly scale with 
xx = 2px/\/~S, i.e., subprocess fractions are similar for 
the same value of xx- 
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FIG. 5: Pseudorapidity dependence of the neutral pion pro- 
duction cross section at = 7TeV for various fixed values 
of pt- The typical theoretical scale uncertainty is indicated 
by the shaded bands for pr = 10 and 50 GeV. 



From panels (a) and (b) one can infer that gluon-gluon 
scattering is the most important channel at relatively 
small pt and central rapidities rj. For pt > 40 GeV at 
i] ~ or at forward rapidities, quark-gluon scattering 
becomes the dominant subprocess. For instance, at large 
forward rapidities, 77 3> 0, one is mainly sensitive to the 
scattering of a valence quark carrying a large momen- 
tum fraction x a off a gluon with Xb <C 1 for which both 
PDFs are large; cf. also panel (f). In the shown kine- 
matic range, quark-quark initiated processes are always 
small and reach a level of 15% only for large px and/or 
large pseudorapidities 77. 

Figures [6] (c) and (d) give the relative contributions of 
quark and gluon fragmentation to the px and 77 differen- 
tial cross sections, respectively. At mid rapidity, gluon- 
to-pion fragmentation is dominant in the entire range of 
Pt shown. For I77I ^> and pt = 25 GeV, quark fragmen- 
tation reaches a level of 80%. We note that for smaller- 
values of pt, the share between quark and gluon induced 
hadron production is roughly equal. 

In the lower panels of Fig. [6] we present estimates of 
the mean momentum fractions {x a ,b) and (z) which are 
predominantly probed in single-inclusive hadron produc- 
tion at the LHC. There are several ways to estimate, for 
instance, an average (z). We define it in the standard 
way by evaluating the convolutions in Eq. ([T]) with an 
extra factor of z in the integrand, divided by the cross 
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FIG. 6: (a) and (b): relative contributions of NLO partonic 
subprocesses da a b initiated by gluon-gluon, quark-gluon, and 
quark-quark scattering to the 7r° cross section at = 7TeV 
as a function of pr and 77, respectively, (c) and (d): same 
as in (a) and (b) but now for the relative fractions of quarks 
and gluons fragmenting into the observed pion. (e) and (f): 
mean values of the momentum fractions x a ,b and z probed in 
7T° production. 



section itself [32|], i.e., schematically we use 



(z) = 



fdzz-^L 

J dzapT 
J dzdpT 



(4) 



(x a ,b) are estimated accordingly. The most important 
thing to notice is that despite the large c.m.s. energies 
available at the LHC, the bulk of the hadrons is pro- 
duced with (z) > 0.4 and (z) > 0.6 for r\ = and 77 3> 1, 
respectively, which is comparable to what one finds at 
much lower c.m.s. energies, for instance, at RHIC. This 
important finding ensures the applicability of the concept 
of FFs which, as discussed in the Introduction, requires 
z > 0.05. Another interesting result is that at mid ra- 
pidity the average x a ,b is considerably larger than what 
one might naively expect from the lower kinematic limit 
x a} b — xt — 2pt/V~S which is of the order a few times 
IO- 3 for pt < 20 GeV. 

Figures [3115] clearly demonstrate the potential impact 
of single-inclusive hadron measurements in a wide kine- 
matic range in pt and rj on future FF global analyses. 
Data will impose stringent constraints on the FFs for 
various hadron species in the large z range at unprece- 
dented large scales \if ~ pt- Scale and to some extent 
also PDF uncertainties are considerably smaller than for 
corresponding measurements at RHIC performed at con- 
siderably smaller values of pt and \/~S which are currently 
used in global fits 0,11]. 

Finally, we discuss the property of xt scaling at LHC 
energies. In pQCD one expects deviations from the naive 
power-law scaling of the invariant cross section a lm in 
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FIG. 7: Illustration of Xt scaling for tv + production at n — 
and LHC energies, (a) and (b): ratios of the invariant cross 
section in Eq. © scaled by p£ for 7 TeV (solid line), 2.36 TeV 
(dashed line), and 0.9 TeV (dot-dashed line) to the result ob- 
tained for y/S = 14 TeV for n = 4 and n — 5, respectively. 
The bands indicate the theoretical uncertainty from varying 
the scales fifj\ r independently in the numerator and denom- 
inator of the ratios, (c): xt dependence of the scaling expo- 
nent n and scale uncertainty (shaded bands) estimated based 
on Eq. ([SJ. 



Eq. © leading to n = n(x T ,VS) d- Figures [7] 
(a) and (b) show ratios of er inv scaled by for 7 TeV, 
2.36 TeV, and 0.9 TeV to the result obtained for y/S = 
14 TeV for fixed n = 4 and n = 5 in Eq. ([2]), respectively. 
As can be seen, scaling violations are sizable for n = 4 
as should be expected from the running of a s and the 
scale evolution of PDFs and FFs. Note that the ratio 
for 0.9 TeV is too large to be displayed in Fig. [7] (a). 
The choice n = 5 in Eq. (J2j leads to much more similar 
ratios for all three energies in the entire range of xt, in 
particular for y/S = 2.36 TeV and 0.9 TeV. 

Power-law scaling of ai nv with a universal exponent n 
can, however, never be more than a rough approximation 
because n must depend on both px and \/~S rather than 
being constant. This kinematic dependence of n is illus- 
trated in Fig. [7] (c) for two different values of px- As in 
Ref. we estimate the scaling exponent n in Eq. ([2]) 
by comparing xt spectra at different c.m.s. energies 
and y/S 7 for fixed pr, i.e., 

\n[a inv (S,x T )/a inv (S',x T )} 

n{XT) = Hv^m — • (5) 



A dependence on both xt and pr is clearly visible and 
needs to be taken into account when comparing theoret- 
ical expectations for xt scaling to experimental spectra. 
Consequently, there is also some ambiguity in estimating 
ti(xt) as the result based on Eq. ([5]) depends also on the 
choice of S and <S". 

In addition, studies of the scaling behavior also suf- 
fer, of course, from theoretical scale ambiguities. The 
shaded bands in Figs. [7] (a) and (b) indicate the uncer- 
tainty from varying the scales Pfj', r by the usual factor 
of two up and down the default choice pr- Compared to 
Ref. we find a much larger variation because, as a 
more conservative choice, we allow for different scales in 
the numerator and denominator of the ratios of invariant 
cross sections shown in panels (a) and (b). 

The possibility of having different factorization scales 
for observables calculated at different c.m.s. energies is, 
e.g. , natural in certain proposed "scale fixing" procedures 
[331 ] . as the "optimum scale" depends on the kinematics 
of the process. It is also interesting to recall recent mea- 
surements of single-inclusive neutral pion production at 
RHIC for three different c.m.s. energies [34{. Although 
the cross sections are described well by pQCD within the 
sizable scale uncertainties, it can be argued that an "op- 
timized" choice of scales in a NLO calculation [ljj would 
suggest to use scales closer to pt/2, pr, and 2p T for 
VS = 62 A, 200, and 500 GeV, respectively. The smaller 
factorization scale at lower c.m.s. energies allows for more 
QCD radiation in the hard scattering matrix elements 
which increases the cross section and leads to a more 
favorable description of the data. To some extent this 
mimics all order resummations of logarithmic contribu- 
tions in the partonic subprocess cross sections which are 
enhanced near the partonic threshold and hence more 
relevant at lower c.m.s. energies for a given value of pt 

B. Unidentified charged hadron spectra 

In the following, we provide a detailed assessment of 
how the present limitations of our knowledge of FFs prop- 
agate to uncertainties for cross section estimates. We fo- 
cus the discussions on unidentified charged hadron spec- 
tra, which, in addition to be accessible with all LHC 
detectors as well as the Tevatron experiments, allow us 
to compare the relative contributions of different hadron 
species and their respective uncertainties. 

Pions clearly dominate the charged hadron spectra, 
with typical uncertainties for the relevant FFs es- 
timated to be at the few percent level j7|- Nevertheless, 
charged kaons, protons, and anti-protons represent non- 
negligible contributions and, potentiated by the much 
larger uncertainties inherent to their FFs 0, H|, may 
yield a significant contribution to the theoretical error 
for estimates of charged particle spectra, perhaps even 
comparable in size to the one for pions. Contributions 
from "residual charged hadrons" , i.e., hadrons other than 
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FIG. 8: Left panel: NLO invariant cross section for the pro- 
duction of charged hadrons for r\ = at four different c.m.s. 
energies. Upper right panel: relative theoretical ambiguity of 
the charged hadron yield at yS = 7TeV induced by varia- 
tions of the scales [ifj> >r an d the uncertainties of the FFs for 
the individual hadron species. Lower right panel: Partial con- 
tribution of the different hadronic species to the the invariant 
charged hadron cross section at y/S = 7TeV. 



pions, kaons, and protons, are known to be completely 
marginal at low energy scales but increase with energy. 
Since the corresponding FFs are very poorly known |8|, 
they can produce sizable uncertainties, comparable or 
even larger than those coming from the much more co- 
piously produced hadron species. Altogether, the com- 
bined theoretical error from FFs may become compara- 
ble to scale uncertainties in certain kinematic regions, in 
particular, at larger transverse momenta. 

To a first approximation, uncertainties from FFs are of- 
ten estimated by simply comparing the results obtained 
with two different optimum fits of FFs. Such a proce- 
dure can give, however, at best a lower bound on the 
true error. In order to arrive at a faithful estimate of 
uncertainties derived from those inherent to FFs, we use 
the robust Lagrange multiplier technique @, which ex- 
plores for any desired observable depending on FFs its 
full range of variations within chosen the range of A\ 2 
tolerated in the fit. Although the actual error estimate is 
more involved than in the standard Hessian method Q, 
it has the advantage of not making any assumptions on 
the shape of the \ 2 profile near the minimum of the fit 
or on how the errors of the fit parameters describing the 
FFs propagate to a given observable. 

As an example, the upper right panels of Figs. [8] and 
[5] show the outcome of propagating the uncertainties of 
the FFs for different hadron species with the Lagrange 
multiplier method to the invariant charged hadron cross 
sections (given in the left panels) for both the LHC and 
the Tevatron at c.m.s. energies of 7 and 1.96 TeV, re- 
spectively. As in the DSS analysis @, Q , we allow for a 
tolerance of A^ 2 /x 2 = 2% in each error analysis. The 



FIG. 9: The same as Fig. [8] but now for pp collisions at 
1.96 TeV. The left panels show also recent data from CDF 
fla l and the ratio of data and a NLO calculation. 



uncertainties resulting from FFs are compared to the the- 
oretical ambiguity from variations of the scales within the 
typical range pr/2 < Hfj>, r < %Pt- Although the scale 
uncertainties are dominant, errors propagated from FFs 
to the invariant charged hadron cross sections are quite 
sizable, and their role becomes increasingly significant at 
larger values of px ■ This is similar to what was found in 
Fig. [H for PDF uncertainties. 

In the lower right panels of Figs. [8] and |9] we show the 
relative contributions of the different hadronic species to 
the inclusive charged hadron spectrum as a function of 
Pt- It is worth noticing that the mixture of hadrons 
is fairly independent of pt, which is linked to constant 
average momentum fraction (z) observed in Fig. |6] (e). 
However, the contribution from residual charged hadrons 
becomes increasingly relevant at larger values of pr but 
still remains fairly small as compared to pions, kaons, 
and protons. 

In Fig. [5] we compare the estimates for invariant cross 
section at NLO accuracy to recent data from CDF [Lsl ]. 
In the lower left panel we also show the ratio of data 
and theory. As can be seen, the theoretical results sig- 
nificantly undershoot the data for pr > 25 GeV which 
caused quite a stir fl9| . Since the data are also in excess 
of the corresponding jet measurement, the most likely ex- 
planation for the observed discrepancy is an experimental 
problem. Nevertheless, it is interesting to know that the 
FF uncertainties estimated above are by far too small to 
account for the hadron yield observed by CDF. 

Finally, it is worth mentioning that when going from 
Tevatron to LHC kinematics there is quite some reduc- 
tion in the relative importance of the residual factoriza- 
tion and renormalization scale dependence as compared 
to FF uncertainties. This again illustrates the impact of 
upcoming LHC data on future global analyses. 
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IV. PION PRODUCTION IN pPb COLLISIONS 

In this Section we propose a set of hadron produc- 
tion measurements in pPb collisions at the LHC as a tool 
to characterize and quantify nuclear modifications in the 
fragmentation process. The LHC is capable of provid- 
ing pPb collisions at a maximum c.m.s. energy of about 
8.8 TeV [27| although such a program is not envisioned 
in the initial phase. However, measurements in dAu col- 
lisions at RHIC were instrumental in interpreting results 
obtained in heavy ion collisions, which are considerably 
more complicated to understand theoretically. 

We note that leading order estimates of neutral pion 
production in pPb collisions at the LHC were recently 
presented in Ref. [36j j but based on the assumption that 
all nuclear modifications can be entirely absorbed into 
nPDFs. In view of the known sizable medium induced ef- 
fects on hadron production yields in, e.g., lepton- nucleus 
collisions [20| , which cannot be explained by nPDFs (2|| , 
such an approach is questionable and may not be ad- 
equate. In fact, in a recent paper j25| the concept of 
medium modified fragmentation functions (nFFs) was 
introduced within the standard factorized framework of 
pQCD. These novel nFFs, which obey ordinary timelike 
scale evolution [r| and have been extracted in a global 
QCD analysis for pions and kaons, allow one to treat hard 
reactions with identified hadrons consistently at NLO ac- 
curacy when combined with nPDFs. This is the approach 
we pursue here to compute predictions for pPb collisions 
at the LHC. 

Conventional factorization of short and long distance 
physics effects is, however, not expected to hold in gen- 
eral in a nuclear environment (Til l2o| , and indeed colli- 
sions between two heavy nuclei show a very large break- 
ing pattern. However, factorization was shown to be 
phenomenologically very successful in describing current 
lepton-nuclei and deuteron-nuclei collision data, includ- 
ing their A dependence (25j . This supports the concept 
of universal, medium modified PDF and FFs at least at 
an approximate level. pPb collisions at the LHC will cer- 
tainly explore the limits of characterizing nuclear modi- 
fications in a factorized pQCD approach. 

To recall the pattern and magnitude of observed 
medium induced modifications for quarks and gluons, we 
show in Fig. [TU] the ratios of nPDFs and nFFs 



Rf(x,Q 2 ) 



ff{x,Q 2 )' t/A 



U ilA 



to the standard PDFs and vacuum pion FFs, respec- 
tively. The ratios are evaluated for Pb at a scale of 
Q 2 = 100 GeV 2 relevant for pion production at trans- 
verse momentum px — lOGeV. 

The l.h.s. of Fig.lTUlshows the nuclear modification fac- 
tors Rf b for u valence quarks, u sea quarks, and the gluon 
g for two standard NLO sets of nPDFs: nDS [22| (solid 
lines) and EPS09 [24| (dashed lines). The deviations of 
Rf b from unity depend on the x region and are typically 




Q" = 100 GeV" 




■ nDS 
EPS09 




10 10 " 10 " 10 



0.2 0.4 0.6 0.8 



FIG. 10: Nuclear modification factors in lead for PDFs [22|,|24( 
(left panels) and FFs [25l ] (right panel) of different flavors at 
Q 2 = 100 GeV 2 as determined in recent global QCD analyses 
at NLO accuracy. 



referred to as "shadowing" (x < 0.01), "anti-shadowing" 
(0-01 < x < 0.2), "EMC effect" (0.2 < x < 0.7), and 
"Fermi motion" (x > 0.7) in the terminology of nuclear 
DIS. Differences between the EPS09 and nDS sets are 
most pronounced for the gluon nPDF, which plays a ma- 
jor role in pPb collisions. Contrary to PDFs, the behavior 
of nPDFs is basically unconstrained below x ~ 0.01 by 
present data, an x region particularly prone to possible 
novel, non-linear features of QCD scale evolution. 

The corresponding medium modifications for FFs are 
displayed in the right panel of Fig. QJJ and are distinctly 
different for quarks and gluons, where one finds suppres- 
sion and enhancement, respectively, compared to vac- 
uum fragmentation functions. The pattern is readily ex- 
plained by the dominant role of quark fragmentation in 
describing the observed hadron attenuation in SIDIS off a 
heavy nucleus, while the enhancement of hadrons in dAu 
collisions is closely linked with the gluon nFF (25|. In 
general, the observed nuclear effects are more pronounced 
for nFFs than for nPDFs, and nFFs can either enhance 
or overturn medium modifications computed with nPDFs 
but vacuum FF as was done, e.g., in Ref. [36| . 

Assuming that factorization holds in pPb collisions at 
a c.m.s. energy of y/S = 8.8 TeV, we present predictions 
for the production of neutral pions in Fig. [TTJ at NLO 
accuracy using nPDFs and the recently proposed set of 
nFFs [25| in Eq. ((JJ. By construction, the entire nuclear 
dependence resides in the non-perturbative nPDFs and 
nFFs, and both the scale evolution and the partonic hard 
scattering cross sections da a b-> c x are taken to be the 
same as in pp collisions. 

The upper left panel shows the pt dependence at cen- 
tral pseudorapidities \rj\ < 0.5 using [if = [ifi = [i r = pr 
in Eq. ((JJ). We refrain from showing alternative results 
for different choices of M/,/',r as the theoretical scale am- 
biguities are very similar to the ones for pp collisions 
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FIG. 11: Transverse momentum (left panels) and pseudora- 
pidity (right panels) dependence of neutral pion production 
in pPb collisions at v5 = 8.8 TeV. The upper panels show 
the cross sections at NLO accuracy, and lower panels display 
the ratio to the corresponding cross section in pp collisions. 
The pt dependent yields are integrated in Wj\ < 0.5 and com- 
puted for dilferent combinations of nPDFs [23, [H, nFFs [25| . 
and vacuum FFs [7(. The rj differential results are given for 
various fixed values of pr, and were obtained using the nDS 
set of nPDFs [H] and nFFs of [1|. 



discussed in Sec. Mil Because of the steep fall of the 
cross section with px over several orders of magnitude 
any differences between the NLO calculations based on 
various combinations of nPDFs and standard vacuum or 
medium modified FFs are hard to notice. The different 
patterns of nuclear modification become clearly visible, 
however, when expressed as ratios to the corresponding 
cross section in pp collisions at the same c.m.s. energy, 
which are given in the lower left panel of Fig. [TT1 

The impact of medium induced effects on the 
hadronization process is clearly visible at larger px by 
comparing the results obtained with the nDS set of 
nPDFs [23 along with nFFs [^J (solid line) and with 
DSS vacuum FFs [7[ (dashed line), leading to hadron at- 
tenuation and enhancement, respectively. The uncertain- 
ties in our present knowledge of nPDFs can be inferred 
from comparing the dot-dashed curve, which uses again 
the DSS vacuum FFs but now the EPS09 set of nPDFs 
24] , with the dashed line. The EPS nPDFs result in a 
somewhat more sizable enhancement than the nDS set 
due to the larger Rg b at medium x values, see Fig. HOl 
which are predominantly probed here; cf. Fig. 1131 below. 
Notice that the EPS09 analysis [II] includes dAu -t tPX 
data from RHIC by assuming, however, that the entire 
nuclear dependence resides only in the initial state, i.e., 
in the nPDFs. This is clearly inadequate in view of 
the sizable hadron attenuation observed in semi-inclusive 
lepton- nucleus scattering [2(| • Hence, in what follows we 
will only use the NLO nDS set of nPDFs, which was also 
obtained in a convolutional approach, very similar to the 



nFF analysis [HI based on Eq. ©. 

The panels on the r.h.s. of Fig. ITT1 show the pseudora- 
pidity dependence of the tt° yields at NLO accuracy for 
various fixed values of px- All results are obtained with 
the nDS set of nPDFs and nFFs of [25]. As before, we 
choose (if = fif = fx r = px in Eq. ([T]). The upper panel 
displays the differential cross sections, and the lower one 
the ratios to the corresponding results in pp collisions. 
The rapidity dependence of the pPb nuclear modification 
factor is of great phenomenological relevance as it probes 
different ranges of momentum fraction Xb in the nPDFs 
while the average value of z in the fragmentation process 
only slowly increases with larger \rj\; cf. Fig. [HI In addi- 
tion, larger values of px at fixed 77 probe both on average 
larger x a ^ and z values. Pions produced in the back- 
ward [r\ < 0) direction of the proton beam require large 
momentum fractions Xb, and nPDFs are mainly probed 
in the ant i- shadowing and EMC region. In contrast, for- 
ward hadrons are sensitive to smaller values of Xf, where 
shadowing is expected; see Fig. [TUJ In general, this leads 
to nuclear modification factors which are asymmetric in 
rapidity for any given value of px ■ 

Based on these kinematic considerations, it was argued 
in Ref. [36| that a detailed study of the px and 77 depen- 
dent hadron yields at the LHC will help to determine 
nPDFs more precisely. However, the neglected medium 
modifications on the FFs can change the theoretical ex- 
pectations for the ratio da p pb/da pp significantly, as can 
be seen in Fig. 1111 In particular, at larger px and rapidi- 
ties our expectations based on consistently including also 
nFFs in Eq. {T]) show much larger nuclear modification 
factors than the few percent effects estimated in Fig. 2 of 
[2(3]. This is readily explained by the strong depletion of 
the quark FFs at large values of z as displayed in Fig. [TU1 
which is most relevant at large 77 and px- 

To elucidate this further, Fig. H"2l gives more details on 
the relevance of different partonic subprocesses (a), the 
role of quarks and gluons in the hadronization (b), and 
the impact of nFFs (c) and nPDFs (d) on the rapidity de- 
pendence of the nuclear modification factor da p pb/da pp . 

Panel (a) illustrates that in the entire px range shown 
in Fig. llll gluon initiated subprocesses dominate the cross 
section at central rapidities. Most of the pions are pro- 
duced from gluon fragmentation as can been seen in 
Fig.[T2"l(b). Only for p T > 50 GeV quark nFFs contribute 
at the level of 20% or more. For hadrons produced at 
larger 77, the fractions in (a) and (b) shift somewhat in 
favor of quark-gluon scattering and quark fragmentation, 
respectively. 

Figures [T2"lfc) and (d) make explicit that the dominant 
medium effect on da p pb /da pp resides indeed in the final- 
state. Panel (c) shows the ratio where the denominator 
is computed with the nDS set of nPDFs such that the 
entire deviation from unity is due to the nFFs. Likewise, 
in panel (d) we compute da pp in the nuclear modifica- 
tion factor with nFFs rather than vacuum FFs. This 
quantifies medium effects due to the initial-state nPDFs. 
The latter effect is much smaller, and our results agree 
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FIG. 12: (a): relative contributions of NLO partonic subpro- 
cesses da a t initiated by gluon-gluon, quark-gluon, and quark- 
quark scattering to the pPb cross section shown in Fig. QT] 
(b): same as in (a) but now for the relative fractions of quarks 
and gluons fragmenting into the observed pion. (c) and (d) 
impact of medium modified FFs and nPDFs, respectively, on 
the rapidity dependence of the nuclear modification factors 
for different values of pr shown in the lower right panel of 
Fig. El 



with the LO estimates shown in Ref. [36|]. Thus we be- 
lieve that single inclusive hadron production at the LHC 
will provide a decisive test of the proposed factorized 
framework for pPb collisions and the concept of medium 
modified FFs. In particular, the ALICE experiment is 
capable of identifying different hadron species in a wide 
range of rapidity which will greatly facilitate such theo- 
retical studies. To disentangle and characterize nuclear 
effects in PDFs and in FFs precisely, requires to study 
also other hard processes in pPb collisions not affected by 
hadronization like jet production, Drell-Yan, or prompt 
photons. Such measurements will provide an important 
input to future global QCD analyses of nPDFs. 

Finally, Fig. ll3l gives an idea which values of x a j, and z 
are probed on average in the pt differential cross section 
at central rapidities shown in Fig. [IT] As for pp collisions, 
the estimates are based on Eq. (Q} . It turns out that mea- 
surements of pion production up to transverse momenta 
of 100 GeV probe both the nPDFs and the usual PDFs 
in the range of momentum fractions from 0.01 to 0.1. To 
access smaller values in the lead nucleus, where novel, 
non-linear features of QCD scale evolution may become 
relevant, one needs to go to very forward pseudorapidities 
[32j . From the lower panel of Fig. [13] one can infer that 
pPb collisions are mainly sensitive to fairly large values 
of momentum fraction taken by the produced pion, with 
(z) slightly increasing with px- This is similar to what 
was observed in dAu collisions at RHIC energies [25| and 
for pp collisions above in Fig. [S] Despite the large c.m.s. 
energy of ^/S = 8.8 TeV, (n)FFs can be safely applied as 
one is fairly insensitive to the region of small z, where the 
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FIG. 13: Upper panel: mean values of x a and Xb probed in 
the proton and the lead nucleus, respectively, for pion produc- 
tion at central rapidities in pPb collisions at %/5 = 8.8 TeV. 
Lower panel: same as above but now for the mean value of z 
computed with the nDS set of nuclear PDFs and DSS vacuum 
FFs 0] (solid line) and nFFS of Ref. HI (dashed line). 



concept of (n)FFs is bound to fail due to finite hadron 
mass effects, higher twist contributions, and the singu- 
lar behavior of the timelike evolution kernels. We note 
that one samples slightly different values (z) in the nu- 
clear and vacuum FFs as can been seen by comparing the 
dashed and solid lines. 



V. SUMMARY AND CONCLUSIONS 

We have presented a comprehensive analysis of single- 
inclusive hadron production in pp and pPb collisions at 
LHC energies based on QCD factorization. 

It was shown that first results from the LHC experi- 
ments for charged hadron spectra agree well with expec- 
tations based on NLO pQCD calculations using latest 
sets of parton distribution and fragmentation functions. 

Based on this success, we have given detailed predic- 
tions for various kinematic distributions of identified and 
unidentified hadrons in pp collisions. Different sources of 
theoretical uncertainties were discussed and estimated. It 
turned out that the residual factorization and renormal- 
ization dependence of the cross sections at NLO accuracy 
represents the dominant source of uncertainty, in par- 
ticular, at small values of transverse momentum of the 
produced hadron. Approximate xx-scaling with n ~ 5 
can be accommodated within the rather large theoretical 
scale ambiguity. Uncertainties from fragmentation func- 
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tions for different hadron species were propagated to the 
single-inclusive yields based on the robust Lagrange mul- 
tiplier method and found to be sizable, although smaller 
than those associated with the truncation of the pertur- 
bative series at NLO. 

To elucidate the possible impact of upcoming hadron 
production data from the LHC on future global analyses 
of fragmentation functions, we have studied in detail the 
relative contributions of different partonic subprocesses 
and the fractions of quark and gluon fragmentation into 
the observed hadrons as functions of transverse momen- 
tum and pseudorapidity. In addition, we have estimated 
the mean values of momentum fractions both in the par- 
ton densities functions and in the fragmentation process 
which are relevant at LHC energies. It was found that 
like at hadron colliders at lower energies, most of the pro- 
duced hadrons take a rather large fraction of the parent 
parton's momentum which ensures the applicability of 
the concept of factorized fragmentation functions also at 
LHC energies. 

Finally, we have proposed a set of measurements of 
single-inclusive hadron production in proton-lead col- 
lisions to shed light on the so far poorly understood 
hadronization mechanism in a nuclear medium. We 
have presented expectations for pion yields at NLO accu- 



racy based on standard QCD factorization using sets of 
medium modified parton distribution and fragmentation 
functions. 

Single-inclusive hadron production in an unprece- 
dented energy range at the LHC will challenge our cur- 
rent understanding of fragmentation functions and help 
to further constrain them, in particular, at large momen- 
tum fractions. The large range of transverse momenta 
of the produced hadrons will allow for detailed studies 
of the scale evolution for fragmentation functions. Data 
obtained in proton-lead collisions will explore the lim- 
its of characterizing nuclear modifications of hadron pro- 
duction yields in a factorized QCD approach and scruti- 
nize the applicability of the recently proposed concept of 
medium modified fragmentation functions. 
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